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ABSTRACT 

We study the relation between the discrete-time version of the flashing ratchet known as Parrondo's games and 
a compression technique used very recently with thermal ratchets for evaluating the transfer of information - 
negentropy - between the Brownian particle and the source of fluctuations. We present some results concerning 
different versions of Parrondo's games, showing all of them a good qualitative agreement between the gain and 
the inverse of the entropy. 
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1. INTRODUCTION 

In the last years, the microscopic field of Brownian particles has received much attention due to new transport 
phenomena under certain conditions. The main features involve breaking of spatial inversion symmetry and 
a thermal bath that drives the system out of equilibrium^ in order to obtain directed transport. The former 
condition is accomplished through an asymmetric potential, usually a ratchet-like potential. The latter can be 
achieved in different ways, either by a periodic or stochastic forcing - what it is known as pulsating ratchets, 
or by an additive driving force unbiased on average - known as tilting ratchets. The flashing ratchet model 
corresponds to the class of pulsating ratchets and it consists on switching on and off either periodically or 
stochastically a ratchet potential. This model has been used recently for DNA transport^ and separation of 
biological macromolecules.^ 

Recently, Arizmendi et. al^ have quantified the transfer of information - negentropy - between the Brownian 
particle and the noncquilibrium source of fluctuations acting on it. They coded the particle motion into a binary 
string of O's and I's, and then zipped the resulting binary file using the Lempel and Ziv algorithm'^ - available 
through the gzip program. In this way they obtained an estimation of the entropy per character /i, as the ratio 
between the length of the zipped file and the original length file, when the text length tends to infinity. They 
used this method to estimate the entropy per character of the ergodic source for different values of the flipping 
rate, with the result that there exists a close relation between the current in the ratchet and the net transfer of 
information in the system. 

The aim of this paper is to apply this technique to a discrete-time and space version of the Brownian ratchet, 
known in the literature as Parrondo's paradox. 

Parrondo's paradox^"^^ is based on a combination of two negatively biased games - losing games - which 
when combined give rise to a positively biased game - winning game. This paradox appeared as a translation 
of the physical model of the Brownian ratchet into game-theoretic terms. But although it was inspired by the 
flashing ratchet model, there has been no quantitative demonstration of their relation until very recently. -'^^"^'^ 

More precisely, the so called Parrondo's paradox is based on the combination of two games (Parrondo's 
games). One of them, say game A, is a simple coin tossing game where the player has a probability p of winning 
one unit of capital and a probability 1 — p of losing one unit of capital. The second game, game B, is a capital 
dependent game, where the probability of winning depends on the capital of the player modulo a given number 
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M . For the original games, M is set to three and there are only two possible values for the winning probabilities, 
Po when the capital is a multiple of three and p\ otherwise. Their numerical values are 



{ Pi= I - e, 

here e is nothing but a biasing parameter that converts games A and B into losing games. When e = it can be 
demonstrated that the fairness condition is fulfilled for both games, that is H^o^Pj = 11^0^(1 ~Pi)- As soon 

as e > the latter condition no longer applies and A and B are both losing games. 

The combination game, game AB, is obtained mixing game A and game B with probability 7. As a result of 
this combination we obtain a winning game, even when e is equal to zero. The corresponding winning probabilities 
are given by the following expressions 

\ P'o= 7P + (1 - 7)Po, /2) 
1 P'\ = 7P+ (1 -7)Pi, 

Several other versions of the games have been introduced. In the so-called cooperative games, ^^'^^ one 
considers an ensemble of interacting players. In the history dependent games, the probabilities of winning 
depend on the history of previous results of winnings and losing. Finally, in the games with self-transition, "'^^ 
there are non null probabilities that the capital can stay the same (not winning or losing) in a given tossing of 
the coins. 

Some previous works in the literature have related Parrondo's games and information theory. Pearce, in Ref.,^'' 
considers the relation between the entropy and the fairness of the games, and the region of the parameter space 
where the entropy of game A is greater than that of B and AB. Harmer et. al^^ studied the relation between 
the fairness of games A and B and the entropy rates considering two approaches. The first one consisted on 
calculating the entropy rates not taking into account the correlations present on game B, finding a good agreement 
between the region of maximum entropy rates and the region of fairness. In the second approach they introduced 
these correlations obtaining lower entropy rates and no significant relation between fairness and entropy rates 
for game B. 

In this paper we aim to relate the current or gain in Parrondo's games with the variation of information 
entropy of the binary file generated. Following the same approach as in Ref.,* we have carried out simulations 
for different versions of Parrondo's games, including the original version, the history dependent games, Parrondo's 
games with self-transition and cooperative games mentioned above. 



2. SIMULATION RESULTS 

We have carried out nmnerical simulations of the different versions of the games. In the simulations, the evolution 
of the capital of the player has been converted to a binary string of O's and I's, corresponding to a decrease and 
an increase of capital respectively. We must remark that each one of these characters is stored as a byte (eight 
bits). A binary file has then been created for different values of the 7 parameter. In order to obtain the desired 
effect there must be a certain delay time, which we will denote by St, between measurements of the capital of 
the player. 

The binary files created have been zipped using the gzip (v. 1.3) program, to obtain an estimation of the 

entropy per character h. For a better comparison with the current or gain - generated when alternating between 
games A and B, we will define R = j^, accounting, in some sense, for the known information about the system. 
In the Appendix there is an extract of the C program used for the case of the original Parrondo's games. 

In Fig.l we compare the gain when alternating between games A and B with probability 7 with the value of 
R obtained for different values of the 6t parameter. The numerical curves are shifted to the origin, that is, we 
plot Ai? = R^ — R-y=0- 

We find a good qualitative agreement between the increase in the gain and the decrease (increase) in entropy 
(i?) as the 7 parameter is varied. Although both curves do not fully coincide, the important fact to consider here 





Figure 2. Comparison between tlie inverse of tiie entropy per ciiaracter R, calculated for different St values, and the gain 
for the Parrondo's games A and B with self-transition varying 7. The values for the probabilities are: p — r = 
po = = 5, Pi = f and ri — j^. See Ref. for an explanation of these parameters. The values used are ot — 6000 

(triangles), 7000 (diamonds), 9000 (squares) and 10000 (circles). 




Figure 3. Curve for the inverse of tlie entropy per ciiaracter R and the gain for the history dependent games, alternating 
between two games with probabilities: Pi = P2 = Ps = |, P4 = j^; Qi = §, 92 = 93 = § and 54 = |. See Ref.^^ for 
an explanation of these parameters. The values for the St are: 3000 (circles), 5000 (squares), 6000 (diamonds) and 7000 
(triangles). 



is that the 7 value for which there is the maximum decrease in entropy agrees with the value for the maximum 
gain in the games. This decrease in the entropy of the system implies that there exists an increase in the amount 
of information that we know about the system. 

In Fig. 2 we compare the gain curve and the entropy curve for the case of Parrondo's games with self- 
transition. Again in this case the maximum gain coincides with the 7 value for the minimum entropy per 
character for all values of St- 

Finally, in Figs. 3 and 4 we show the comparison in the case of the history dependent games, and cooperative 
games, sharing all them the same features as the previous games analyzed before. 

3. THEORETICAL ANALYSIS 

In this section we develop a simple argument that explains the observed relation between the gain of the games 
and their entropy. In the seminal work by Shannon^^ the entropy per character of a text produced by an ergodic 
source is given by the following expression 

F=-^p,-l0g(p,;) (3) 

i 

where pi denotes the probability that the source will 'emit' a given symbol a^, and the sum is taken over all 
possible symbols that the source can emit. 

If the source can be found in several states, each one with probability Pj, then the entropy reads 

M 

H = Y,P,H, = -Y,P,pilog(j4). (4) 

j = l j,i 

In the last equation pj denotes the probability of emitting the symbol when the source is found in the state j. 
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Figure 4. Curve for the inverse of tlie entropy per ciiaracter R and tiie gain for the cooperative version with probabilities: 
p = ^, pi = 1, P2 = P3 = P4 = and A'' — 150 players for different values of the 7 parameter. See Ref.^^ for an 
explanation of these parameters. The values of the St are: 6000 (circles), 8000 (squares), 12QQQ (diamonds) and 18000 
(triangles). 



Remember that we measure and store the capital of the player after St rounds have been played. When 6t 
equals one it means that we measure the capital of the player in each round. Translating this into entropic 
terms we have that, considering the original Parrondo's games, our source will be in two possible states, say, 
when the capital is multiple of three and otherwise. From Markov chain theory^^ we know that the stationary 
probabilities of finding the capital in each of the states are given by 

Uo = j^il-pi+pl) 

(5) 

Hi = l-Ho = ^{2-po-pi + 2poPi) 
where D ^ 3 - po - 2pi + 2poJ3i + pI- 

It is interesting to consider our source (the capital of the player) as composed of several sources, each one 
with an homogeneous statistical structure, i.e. they are ergodic sources. Each source then will have a determined 
probability of 'emitting' a zero or a one, as the player has a different winning probability depending on whether 
his capital is multiple of three or not. If we consider that the source will emit the symbol 1 with probability pi 
and with probability 1 — Pi then Eq. 4 reads 



H =i:^i^,H, = -j:^^n,p>\og{pi) 

= -Ho ■ {pi ■ log(pi) + (1 -pi) ■ log(l -pi)) - (1 - Ho) • (P2 ■log(p2) + (1 -P2) ■ log(l -P2)) 

In Fig. 5 we plot the inverse of Eq.6 for the values of p^ and 11^ corresponding to the original Parrondo's games 
together with the numerical curve obtained when zipping the files obtained with 6t = 1- It is worth noting how 
the negentropy decreases as 7 increases, translating into an increase of known information about the system. 

Alternatively, we can now estimate how the entropy changes for (5f 3> 1, when we allow a large amount of time 
between measures. As St is increased, the statistical description used above cannot be used anymore. For higher 
St values, the system gradually loses its memory about its previous state and so what we obtain are measures 




Figure 5. Plot of the theoretical curve for the inverse of the entropy per character and the numerical curve (circles) for 
the case when St = 1 step. 



that can be considered as statistically independent, as we will not know for sure which will be the probability 
of finding the source in any of the two possible states. Now the dominant term will be the averaged probability 
-over all possible states- of finding the capital of the player with a value higher than the previous measure. 

Considering the same stationary probabilities Hq and Hi given by Eq.5 we can obtain the expressions for the 
average winning and losing probabilities over these two states as 

Pwin = Y^jii nj Pj =UoPo + Hi p[ 

It can be noticed in all Figs. 1,2, 3, 4 that as St increases, the entropy curve acquires a more pronounced peak. 
This effect can be explained due to the fact that as the time between measures increases, the probability of 
finding the player with a higher (lower) capital increases (decreases). We can give an estimate of the probability 
that after St rounds played, the capital of the player will be greater (p>) than it was initially by considering the 
probability that there have been less that St/2 loses in the St coin tosses, 
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These previous equations together with 

if = •log(pj-(l-p>)-log(l-pj, (9) 

allows one to obtain the entropy per character. In Fig. 6 we can see the curve for the variation of the inverse of 
information entropy for different values St- In this case the negentropy decreases as St increases, reflecting an 
increase of known information about the system. 

This theoretical explanation also holds for the case of history dependent Parrondo's games and Parrondo's 
games with self-transition. For the latter case, though, we must take into account the fact that there exists a 
non negligible possibility that the capital can remain the same after a round played. In this case we need only 
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Figure 6. Plot of the theoretical curves for the inverse of the entropy per character for different values of St, together 
with the numerical curve obtained for St = 500 steps. 

to consider a source that will emit three different symbols, depending on whether the capital is greater, equal or 
smaller than our previous measure. 

As a conclusion, we have demonstrated that there exists a transfer of information in the case of Parrondo's 
games, considered as a discrete-time and space version of the flashing ratchet. This effect takes place in every 
existing version of the games analyzed, showing the robustness of it. The quahtative analysis for the case of the 
original Parrondo's games have demonstrated the fact that the capital of the player can be thought of as a mixed 
source, composed of two ergodic sources, depending on whether the capital is a multiple of three or not. 

We have also given some insight about the behavior of the information entropy for low and high values of 5t , 
and why there is a decrease (increase) in the entropy per character h (information) as 5t increases. 
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APPENDIX 

We provide the C program used for creating the binary files that will be zipped for obtaining the entropy per 
character h. 

int mainO 
{ 

int i, j , k, ijk; 
int ijkl , i j , n; 

long lunax, irepet, icont, deltat, t; 
long Nrepet , c, cO; 
double gcunma, pa, pbl, pb2; 
double u, V, gain; 
unsigned char cero=0; 
unsigned char uno=255; 



srand48 (342342) ; 



FILE *current; 

current = f open ("gain. dat" , "w") ; 
/* We fix the lenght in bytes of the binary file */ 

Nrepet=300000; 
/* We impose the delay time between measurements */ 

deltat = 500; 



/* Probabilities of the original Parrondo's games */ 

pa = 1./2. ; 
pbl = I./IO.; 
pb2 = 3./4. ; 



ijk = 0; 

for (gamma = 0.; gamma <= 1.0; gamma = gamma + 0.04)-[ 

/* we create a binary file for each value of gamma */ 

sprint f (q, "gamma7od.dat" ,ijk) ; 

FILE *data; 

data = f open(q, "w") ; 

++ijk; 

cO = 0; 
c = 0; 
i= 0; 

ijkl = ; 
t = 0; 

while (i<Nrepet){ 
++t; 

if (draiid48() < gamma ){ 

/* we play game A */ 

if (drand48() < pa ){ 

++c;} 

else { 

~c; 

} 

++ijkl;} 
else{ 



/* we play geime B */ 

if (fmod(c,3) == ){ 

if (drand48() < pbl ){ 

++c;> 

else { 

--c;} 

} 

else { 

if (drand48() < pb2 ){ 

++c;> 

else { 

--c;} 

} 

++ijkl;} 

/* every time that ijkl equals deltat we evaluate the actual capital 
of the player compared to its previous value */ 

if (ijkl == deltat) {. 
ijkl = 0; 

if ((c-cO) > 0) { 

putc(cero,data) ;> 

else{ 

putc(uno,data) ; 
} 

i++; 
cO = c; 
} 
} 

/* we also write the gain for every gamma value */ 

gain = (double) c/( (double) t) ; 

f printf (current , "°/.f\t°/,f \n" , gamma, gain) ; 

f close (data) ; 

> 

f close (current) ; 
return EXIT_SUCCESS ; 
> 
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